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Abstract 



Radiative corrections to the Zemach contribution of the hydro- 
gen hyperfine splitting are calculated. Their contributions amount 
to — 0.63(3) ppm to the HFS. The radiative recoil corrections are es- 
timated to be 0.09(3) ppm and heavy particle vacuum polarization 
shifts the HFS by 0.10(2) ppm. The status of the nuclear-dependent 
contributions are considered. From the comparison of theory and ex- 
periment the proton polarizability contribution of 3.5(9) ppm is found. 
The nuclear structure-dependent corrections to the difference f/i/ s (ls)— 
n 3 Vhf s (ns) are also obtained. 
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1 Introduction 



The hyperfine splitting of the ground state of the hydrogen atom is one of 
the most precise measured values fj], || 



uhfsO-s) = 1420405.7517667(9) kHz, (1) 

but the theory is not able to obtain this result with such an accuracy The 
main problem is due to the proton structure. This work is devoted to nuclear 
structure-dependent contributions. 

An expression of the leading nuclear structure-dependent correction to 
the hyperfine splitting in the hydrogen atom was found by Zemach || . Later 
nuclear structure-dependent corrections were investigated in Refs. ||, ||, H, 
|7|, [|, but the radiative structure-dependent corrections have not yet been 
found. The pure radiative corrections are evaluated in the next section. 
The third section is devoted to radiative recoil corrections. After that the 
vacuum polarization of heavy particle is discussed. The last sections consider 
the status of the nuclear structure-dependent terms. 



2 External Field Approximation 

2.1 Leading structure-dependent term 

The Zemach expression of the leading nuclear structure-dependent correction 
to the ground state hyperfine splitting in the hydrogen atom has the form || 

Au(Zemach) = vp (—2Zam e J J d 3 r d 3 r' pe(t)pm(y')\y — r ' 

or 

2Zam e r d 3 p 



Av(Zemach) = vp 



7T 2 J p 4 



G E (p 2 )G M (p 2 

l + K 



(2) 



where Pe{ y ) and Pm(f) are the proton electric charge and magnetic moment 
distribution respectively G_b/a/(p 2 ) is the Sachs electric/magnetic form fac- 
tor, k is the proton anomalous magnetic moment. The Fermi energy vp is 
defined as the ground state hyperfine splitting in the nonrelativistic theory. 
Here, relativistic units in which % — c — 1 and a = e 2 are used. Z is the 
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nuclear charge in units of the proton charge. It is equal to one in the hydro- 
gen atom, but some results in this paper like eq.@ are valid for any low-Z 
hydrogen-like atom. 



2.2 Dipole approximation 

In the well-known dipole approximation of the electromagnetic form factors 
of the proton 



G E (p 



2, G M ( P 2 ) 



1 + K 



it is easy to solve eq.( 



Au(Zemach) = vp 



A 2 



A 2 + p 2 



35 Zam t 
"T A 



(3) 



(4) 



2.3 Vacuum polarization 

To calculate the contribution of the electronic vacuum polarization one can 
insert the well-known asymptotic behaviour of the polarization operator 



a 



71 



-logP--- 

3 ml 9 



(5) 



into the right hand side of eq.@. Integrating eq.fl2|) within the dipole 
approximation(eq.(|3])) yields: 



Av '(structure — VP) = Au(Zemach) 



a 



71 



2, A 2 
o lo S- 



mi 



634 
315 



This result is in agreement with an estimate found in Ref. 

The logarithmic part of this result may be used for any low-Z hydrogen- 
like atom (if the magnetic square radius (Rm) is aproximately equal to the 
charge square radius (Re)) after substituting A — > \/Yl/R E . 
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2.4 Self energy 



The self-energy contribution is evaluated by using an explicit asymptotic 
expression of the one-loop insertion into the electron line. The result is 

Au( structure — e-line) = AviZemach) ■ — ( — X . (6) 

7T I 4J 

The coefficient (-5/4) arises from nontrivial radiative insertion into the elec- 
tron line (-7/4) |§ and from the anomalous magnetic moment contribution 
(1/2). It should be mentioned that this result has been obtained without the 
use of the dipole fit (eq.(|3|)) and it is also valid for any low-Z hydrogen-like 
atom. It can be used in a wide interval of the nuclear charge up to Z = 25 
and the uncertainty is expected as to grow (irZa) 2 in relative units. 



2.5 Binding corrections 

The higher order binding corrections have the order {Za) 2 or (Za m e )/A in 
units of Av(Zemach). They are small and it is enough to take into account 
only their logarithmic parts. The (Za) 2 -term arises from the Dirac correction 
to the wave function 

(Za) 2 1 
Av(structure — Dirac) = Av(Zeraach) ■ log 



(Za) 2 ' 

The (Zam e )/A-contribution is due to the nuclear charge distribution and 
the Fermi interaction. The result is 

/ 2 1 \ 

Au(structure — charge) = Avp ■ —-(Zam e ) 2 r 2 \og 



(Zaf 



These corrections are small, but they are considered here, because only 
the binding corrections can contribute to the hfs splitting of the higher-/ state 
or to the difference 

A hfs (n) = v hf8 (ls) - n 3 v hfs (ns). 

The result for the difference can be found in the non-relativistic approx- 
imation (cf. Ref. |I4"|) 
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SAhjs^n) = Ais(Zemach)-(Zay 



, s , in - l)(n + 9) 

i/;(n+ 1) --0(2) -logn- - 



4n 2 

(7) 

where VK 2 ) = (d/dz) \ogT(z), and for the hfs of states with / > 
. , . . Av(Zemach) 2 n 2 - 1 

<^K) = ^ y Za ) -^-Wi- 

The (Zo! m e )/A-contribution may be found in the same technics useds in Ref. 



2,_ . 2 2 Wn- 1 



5A h/s (n) = Ai/ F - (^--(Zam e ) 2 r 2 J — - logn + ^(n) - V(1)J (8) 

The binding corrections of eq.(|7]) and eq.(|]) shift the difference A^/ s (2) by 
10~ 9 kHz only. 



3 Recoil Contributions 

3.1 The leading term 

The final result of the pure recoil corrections was found in Ref. pj by nu- 
merical means using the dipole approximation. It includes large numerical 
cancelation between different terms. This cancelation can be understood 
analytically from the leading logarithmic term j|, [5[ for atoms with a non- 
structured nucleus (terms VO, VV and k 2 of Ref. M) 

3Zam m p 2(1 + k) - (1 + k) 2 + \k 2 
Au(rec) = —v F ■ log — - • — . (9) 

7r m p m e 1 + k 

The coefficient 

2(1 + k) - (1 + k) 2 + \k 2 

1 + K 

is equal to 1 in muonium and 0.070... in the hydrogen atom. 

In the following subsections radiative recoil corrections (i. e. radiative 
corrections to eq.(pl)) are considered. 
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3.2 Electronic vacuum polarization 



The electronic vacuum polarization term contains the same structure of the 
Dirac matrix as the leading term of eq.(§) and the same cancelation occurs. 
The estimate 



Av{VP - log) = -v F lo § 



2a(Za) m P , 9 A 
^—\og 2 — 



2(1 + k) 



3,^2 



- (1 + KY + fK 
1 + K 



can be obtained by using the asymptotics of eq.(|5|) (see also 



3.3 Self energy 

The self-energy recoil contribution can be estimated from this contribution in 
the muonium atom, which includes a non-relativistic pole ("the J'-term" in 

10| or "the iVi?-contribution" in Refs. 



11, 12 



the definitions used in Ref. 
logarithmic and constant contributions of two different structures of the Dirac 
matrix [10, |ll], [L2|. The numerically important contributions arise from the 
pole and logarithmic terms. They can easily be adjusted to the hydrogen 
atom. The constant can be used to estimate unsertainty. The anomalous 
magnetic moment contribution has to be added as well. The final numerical 
result is presented in the last section. 



4 Heavy Particle Vacuum Polarization 

The correction due to muonic and hadronic vacuum polarization have been 
treated for a point-like nucleus in Ref. ||16|f |. These corrections are found 
here in the external field approximation (for details see Ref. [16}|) with the 
dipole form factors of eq.(£J). The numerical results are presented in the last 
section. The results for the point-like nucleus 



1 That work contains some misprints: the left part of eq. (9) should be multiplied by 
2ir; the result for the hadronic contribution to the HFS in eq.(10) and in the Table should 
be multipied by 2. 
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3 7Yt 

Au(u — VP — point) = — a(Za) — - vf, 

4 

and 

Au(hadr — VP — point) ~ 0.7(3) Ai/(/x — VP — point) 

and for the finite- size nucleus are different. The finite-size nucleus results are 
only some the 30% of the point-like nuclear corrections. The results for the 
finite-size nucleus are given in this work within the external field approxima- 
tion. The uncertainties are estimated by the unknown recoil contributions 
(cf. Ref. TO). 



5 Parameterization of the Dipole Fit 

Because no reliable self-consistent values of the Zemach correction are known 
numerical results can be obtained only after reconsideration this correction. 
The parameter A which is needed for this calcilation is directly connected 
with the proton square charge radius 
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p 



A 



Comparision of recent experimental results with the theory (see e. g. Ref. 
H]]) favors the newer value JTBj 



r p = 0.862(12) fm, 



or 



A = 0.845(12) m p . (10) 
However, as most of the work in this field was done more than 15 years ago, 



the older proton radius of Ref. [20| 

r p = 0.809(11) fm, 

or 
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A = 0.898(13) m f 



(11) 



was used. It is well-known that the form factor from the older radius and 
the value A = 0.898(13) m p is good as long as the momentum transfer is 
not too low. Hence, the first problem is to understand what momenta are 
important in the integration of eq.(||]). In order to solve this problem high 
and low momenta have been separeted in the integral 

a i\ 8Zam e 

Av(Zemacli) = v F x 

7r m„ 



m. 



Q dp 



- 1 



dp 
Iq p 2 



;i2) 



The low momenta asymptotics has the form 



(G D (p 2 



A 2: 



and the high momenta asymptotic behaviour is 



(13) 



(Gd(p 2 ) 



1 



-1. 



(14) 



The results of integrations are presented in Fig.l. 

The main radius-dependent contribution arises from the low momenta 
asymptotics (— 4p 2 /A 2 ), which should be directly connected to the radius. 
Our results for several Q (see also Fig. 2) 





' -40.92(59) 


10~ 6 


• v F , 


Q 


= 0.30 m p , 


Ai>(Zemach) = < 


-41.07(68) 


10~ 6 


■ VF, 


Q 


= 0.35 m p , 




t -41.24(68) 


10~ 6 


■ v F , 


Q 


= 0.40 mp, 



are obtained as the sum of the 4p 2 /A 2 -contibution from the lower momenta 
with A = 0.845 m p and the average value of the remaining terms with A = 
0.845 m p and A = 0.898 m p . The uncertainty is obtained from the sum of 
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the squares of the parameter-induced uncertainty in the 4p 2 /A 2 term and 
half of the difference for the contribution of G 2 D — 1 — 4p 2 /A 2 and of higher 
momenta with A = 0.845 m p and A = 0.898 m p . The details of calculations 
are contained in the appendix. 

The results for different values of Q are almost independent of Q. They 
should be compared with results for the point-like proton 



Ai>(Zemach) 



' -41.15(58) 


io- 6 


• v F , 


A 


= OA 


345(12) m p , 








r P 


= OA 


362(12) fm 


-38.72(56) 


10~ 6 


■ VF, 


A 


= OA 


398(13) m p , 








r p 


= OA 


309(11) fm 


-40.43(43) 


10~ 6 


■ VF, 


A 


= 0. 


860(9) m p , 








r p 


= 0. 


847(9) fm. 



(16) 



The last proton radius value is the result presented in the recent work |2T 
One can see that the results of our estimate in eq 



Iq) are close to the 
result of eq.fllBD for A = 0.845(12) m p , but the uncertainties are a little 
higher. This agreement of eq.(|15|) with eq. (|T6|) for A = 0.845(12) m p is due 
to the use of this value in our calculation of the 4p 2 /A 2 term. Within the 
interval between 0.3 m p and 0.4 m p of value of Q the approach described here 
is expected to yields the best results. The value 



Av(Zemach) = -41.07(75) • 10" 6 • v F 



(17) 



is used for futher numerical calculation. We expect that this result is more 
safe than any simple dipole fit values from eq.(|IT 



6 Conclusion 

In one of the latest works |§ , devoted to the ground state hyperfine splitting 
in the hydrogen atom, compararison of theory and experiment leads to the 
difference 



VHFsjexp) - v H Fs(theo) 
VHFs(exp) 



(0.56 ± 0.48) ppm. (18) 
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The theoretical exression exludes the unknown proton polarizability so it 
may be estimated by the difference in eq. fll8D . The theoretical limitation for 
the proton polarizability contribution is fl2"2"|, E31 



\5(polarizability)\ < 4ppm. (19) 
The result of this work is 

VHFsjexp) - v HFS {theo) , x 

r = 3.5 ± 0.9 )ppm (20) 

VHFs{exp) v ' 

instead eq. ( |T8| ) . The changes of the theoretical values are presented in Tables 
1-3. The older comparission of theory and experiment in eq.(|18D implies that 
the polarizability contribution is much lower than the limitation of eq. ([19|). 
However, our comparision in eq.(p0|) leads to a result close to this limitation. 



The hyperfme splitting of the hydrogen ground state is more sensitive to 
the proton structure value than the Lamb shift. We hope that investigation 
of the hfs will lead to a better understanding of the proton and to a more 
accurate calculation of the Lamb shift, the precision of which is limitted by 
the proton radius. 
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Planck-Institut fur Quantenoptik and he is very grateful to T. W. Hansch for 
his hospitality. The author would like to thank Thomas Udem for reading 
the manuscript and useful remarks. 
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sian Foundation for Basic Research. The presentation of the work at the 
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A Q-integrals for calculation of the Zemach 
correction 



1. The integrals over different area are 
I < (Q,A)=m p J o Q ^ 



(G D (p 2 ) 



35 m p 
— — — — arctg I — 



16 A 



Q\ 19 m p Q 11 m p QA 2 1 m p QA 4 



and 



A J 16A 2 + Q 2 24(A 2 + Q 2 ) 2 6 (A 2 + Q 2 ) 3 



I > (Q,A)=m p r^\(G D (p 2 ) 



35 m r 



'A' 



— —arctg \ - ] - 



19 m p Q 11 m p QA 2 



1 m p QA 4 



QJ 16A 2 + Q 2 24 (A 2 + Q 2 ) 2 6 (A 2 + Q 



2\ 3 ' 



2. The contributions of asymptotics of eq.([L3|) and eq.([L4|) are easy to find: 

A<(Q,A) = -4 

and 



m p Q 



A > (Q,A) = - T ^. 

The contributions of remaining low-momenta term is denoted by 

R < (Q,A) = I < (Q,A)-A < (Q,A). 

All results with different A are presented in Fig. 1 as func- 
tions of Q. 

3. The following combinations are used as the result and the uncertainty 
I(Q) = A <( Q, 0.845 m p ) + R<(Q-°**>»,H UQASXm,) 
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I>(Q, 0.845 m p ) + I>(Q, 0.898 m p ) 



and 

2 



5I{Q) = 



(85Am p Q^ 
\ A3 t 



+ ^i?<(Q, 0.845 m p ) -R<(Q, 0.898 m p )y + ^7>(Q, 0.845 m p ) - 7>(Q, 0.898 m p )^ 2 ^ 1/2 
The function /(Q) is presented in Fig. 2. 
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List of captions to the tables and the figures 



Table 1: Old and new values of some relative contributions to the theoretical 
hfs (in ppm). 

a The following new values have been used: a -1 = 137.0359994(6) from Ref. 
Q and m p /m e = 1836.152667(4) from Ref. |2 . 

b The result of Refs. |33], [!4j used in Ref. is incorrect (for details see 
and Refs. 



review 



35, 36 



c This is actually a result of Ref. |J, but the older proton radius was used 
there. The uncertainty given there was tripled here. 



Table 2: Contributions to the relative shift of the theoretical result for the 
hfs (in ppm). 

Table 3: Contributions to the uncertainty of the theoretical result for the hfs 
(in ppm). 

Fig 1: Contribution to the integral in eq. ([12]) for A = 0.898 m p . J> - con- 
tribution from higher momenta, A < - contribution of asymptotics of eq.(|T3|) 
from lower momenta and i?< contribution of the remaining terms. 



Fig 2: The Zemach correction as a function of Q. 
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Term 


Old 


New 


Ref. 




-1102.15 


-1102.28 




a 






a(Za) 2 


1.90(4) 


1.84(12) 




b 






a 2 (Za) 




0.09 






, 26, 


27 




higher ord. 




0.01(2) 




28, 


2f 


1 


VP — Structure 




-0.74(1) 


this work 


SE — Structure 




0.12 


this work 


Binding — Structure 




-0.01 


this work 


Zemach 


-38.72(56) 


-41.07(75) 


this work 


Recoil — Structure 


5.22(14) 


5.22(42) 




c 






Recoil — VP — Structure 




-0.02 


this work 


Recoil — SE — Structure 




0.11(2) 


this work 


\i — VP — Structure 




0.07(2) 


this work 


Hadr — VP — Structure 




0.03(1) 


this work 


Weak interaction 




-0.06 




|30| 





Table 1: 
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Term 


Shift 


Constants (a ,m e /m p ) 


-U.lo 


QED 


0.04 


Zemach 


-2.36 


Radiative — Structure 


-0.63 


RRC — Structure 


0.09 


Heavy — VP — Structure 


0.10 


Weak interaction 


-0.06 


Total 


-2.95 



Table 2: 



Term 


Uncertainty 


QED 


0.12 


Zemach 


0.75 


Radiative — Structure 


0.03 


Recoil 


0.42 


RRC — Structure 


0.03 


Heavy — VP — Structure 


0.02 


Total 


0.87 



Table 3: 
19 



I=i<3 1 




Fig. 2 



-38.5 I 1 1 1 1 1 1 1 1 1 1 1 r 



-39.5 - 

o 




-43.5 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.05 0.15 0.25 0.35 0.45 0.55 0.65 



Parameter Q 



